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abstract. Sequence data of the chloroplast matK gene generated a phy- 
logeny of Liquidambar containing two robust clades. One clade consisted of 
the Chinese species L. acalycina and L. formosana, while the other was com¬ 
posed of L. orientalis from Turkey and the North American L. styraciflua. 
The data support a close relationship between the western Asian and North 
American species, but not the division of Liquidambar into section Cathay- 
ambar (L. formosana) and section Euliquidambar (L. acalycina, L. orientalis 
and L. styraciflua). Sequence divergence of the matK gene ranged from 0.1 
to 1.0% among Liquidambar species and the estimated divergence times of 
the disjunct species in the genus were 45-90 mya, which agrees with the 
fossil record. 
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Liquidambar L. is the only genus in the Hamamelidaceae that 
has a disjunct distribution with species occurring in western Asia, 
eastern Asia, and North America. Four species are currently rec¬ 
ognized in the genus. Liquidambar formosana Hance is wide¬ 
spread in eastern Asia (Chang 1979; Li 1977). Liquidambar aca¬ 
lycina Chang is found in at least nine provinces in mainland Chi¬ 
na (Chang 1979). Liquidambar orientalis Mill, occurs in Turkey 
and some nearby islands such as Rhodes (Rechinger 1943) and 
Cyprus (Holmboe 1914). Liquidambar styraciflua L. is widely 
distributed in eastern and southeastern North America and south¬ 
ward at high elevations in the mountains ot Mexico and Central 

America to Honduras (Bogle 1986). 

In Liquidambar, L. acalycina and L. formosana ditter from L. 
orientalis and L. styraciflua in having 3-lobed instead of 5 (-7)- 
lobed leaves. The leaf lobes of L. orientalis are further subdiv id- 
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ed, but they can intergrade with those of L. styracijiuci (Bogle 
1986). The presence of “setae” (Harms 1930), the carpel-like 
organs in pistillate flowers of L. formosana, has been used to 
distinguish L. formosana as section Cathayambar Harms from 
the other Liquidambar species as section Euliquidambar Harms 
(Harms 1930; Chang 1979). 

Disjunct distribution of closely related plant species has long 
attracted the attention of both plant systematists and biogeogra¬ 
phers (see review in Boufford and Spongberg 1983; Crawford 
and Lee 1992; Lee and Crawford 1991; Lee et al. 1996; Tiffney 
1985a, b; Wen et al. 1996; Wen and Zimmer 1996). For Liquid¬ 
ambar, Hoey and Parks (1991, 1994) studied genetic divergence 
ol the four species using allozymes and found that L. orientalis 
was more closely related to L. styraciflua than to either of the 
two eastern Asian species. Their study provided support for the 
existence of Atlantic land bridges between eastern North America 
and western Europe in the upper Cretaceous and Tertiary periods. 

Crawford et al. (1992) pointed out that DNA sequence data 
may provide a more precise estimate of divergence than second¬ 
ary chemicals or allozymes. Recent studies support that sugges¬ 
tion, and have shown that both nuclear and chloroplast DNA 
sequence data are informative in resolving phylogenetic relation¬ 
ships of disjunct taxa, even though resolution is variable for dif¬ 
ferent taxa at different levels (Kim and Jansen 1994; Suh et al. 
1993; Xiang et al. 1994). Undoubtedly, a high resolution of phy¬ 
togeny provides a foundation for understanding phytogeography 
of disjuncts (Wen and Zimmer 1996). 

The objective of this study is to use DNA sequences of the 

chloroplast matK gene to investigate genetic divergence of the 

species of Liquidambar at the nucleotide level and to evaluate 

the phylogenetic and biogeographic relationships among the four 
Liquidambar species. 


MATERIALS AND METHODS 

Fresh leaves were collected from small trees of Liquidambar 
styraciflua, L. acalycina, and L. formosana cultivated in the 
greenhouse of the University of New Hampshire. Leaf buds of 
L. orientalis were provided by Tracy Omar at the University ot 
Washington Arboretum, Seattle. Leaves of Mytilaria laosensis 
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Table 1. Locations and base compositions of amplification and sequenc¬ 
ing primers used in this study. * this primer was synthesized with equal parts 
of “C" and “T” at base position 6. 

Primer 5’ sequence 3’ Designed by 

Forward 


matY>F\ 

matY^Fl 

m«rKF4 

matKFS 

niarKF6 

Reverse 

marKRl 

matKR2 

matKR2-2 

»iatKR3 


ACT 

GTA 

TCG 

CAC 

TAT 

GTA 

TCA 

Tao Sang 

GTT 

CAC 

TAA 

TTG 

TGA 

AAC 

GT 

Tao Sang 

ACC 

CCA 

CCC 

CAT 

CCA 

TCT 


Jianhua Li 

TGG 

AGY 

CCT 

TCT 

TGA 

GCG 

A* 

Jianhua Li 

TCA 

GTG 

GTA 

CGG 

AAT 

CAA 

ATG C 

Jianhua Li 


GAA 

CTA 

GTC 

GGA 

TGG 

AGT 

AG 

Tao Sang 

TTC 

ATG 

ATT 

GGC 

CAG 

ATC 

A 

Tao Sang 

ACG 

GGG 

CCA 

TAA 

GAA 

AGT 

CG 

Jianhua Li 

GAT 

CCG 

CTG 

TGA 

TAA 

TGA 

GA 

Tao Sang 


Lee. were provided by Zhong-chun Luo at the Forest Bureau of 
Xinning, Hunan, China. 

Total genomic DNAs were extracted from fresh leaves or buds 
following the protocol of Doyle and Doyle (1987). Polymerase 
Chain Reaction® (PCR) was conducted in 0.2 ml thin-walled 
microcentrifuge tubes. Each 50 p.1 reaction included 5 pi of 10x 
Taq extender buffer (Stratagene, CA), 4 p.1 of 2.5 pm dNTP, 4 
units (0.8 p.1) of Taq extender (Stratagene, CA) and Taq poly¬ 
merase (Promega, WI), 1 p.1 20 p.m primers, 2-3 p.1 genomic 
DNA solution (50-100 ng DNA), and an appropriate amount of 
UV-treated distilled water. The PCR thermocycler program fol¬ 
lowed Johnson and Soltis (1995) and the primers were matKFl 
and /raj/KRL The PCR products were loaded on 0.8% LMP (Low 
Melting Point) agarose gel along with lambda Hindlll i >NA size 
markers and run for 2—3 hours at 40 volts in 0.5x TBE butter. 
The band identified by comparison to the markers was then ex¬ 
cised from the gel, liquefied at 65°C, and treated with agarase tor 
30 minutes at 37°C. This gel-purified PCR product was used di¬ 
rectly as a sequencing template. 

Sequencing reactions were carried out using the Cycle Se¬ 
quencing Kit and following the manufacturer s protocols (Applied 
Biosystems, CA). The primers for sequencing were wafKFl, 
mat KF2, matKF4, matKl 5, matK.F6, /no/KRl, /?mfKR2, 
wa/KR2-2 and w«/KR3 (Table 1). The approximate locations and 
exact base compositions of the matK primers are shown in Figure 
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matK FI F2 F4 F5 F6 R3 R2 R2-2 

-> -> -> -> -> <- <- <- <- R1 

Figure 1. Approximate relative locations of matK primers (base compo¬ 
sitions are listed in Table 1, shaded areas are introns). 

1 and Table 1. The cycle sequencing products then were separated 
on 6% polyacrylamide gel using an Automated Sequencer 373A 
(Applied Biosystems, CA) in the Sequencing Facility Center at 
the University of New Hampshire. The sequence chromatograms 
were analyzed using the SEQED program (Applied Biosystems, 
CA). Also, in order to assure correct basecalling, we overlapped 
sequences generated from adjacent primers of either the same or 
opposite directions. The complete matK sequences of Saxifraga 
integrifolia (GenBank accession number L20131) and Sullivantia 
siillivantii (GenBank accession number L20130) were used to de¬ 
termine the limits of matK sequences in Liquidambar. 

The sequences were readily aligned using the MEGALJGN 
program of DNA* software packages (DNA* Inc., WI) and by 
sight. The aligned sequences were imported into the PAUF (Phy¬ 
logenetic Analysis Using Parsimony) computer program (Swof- 
ford 1993) to search for the shortest trees, using the exhaustive 
search option. Our analysis of intergeneric relationships of the 
Hamamelidaceae using sequences of internal transcribed spacers 
(ITS) of nrDNA has shown that Mytilaria laosensis is the sister 
taxon to Liquidambar (Li et al., unpubl.). Thus, M. laosensis was 
used as the outgroup in this analysis. All characters and their 
states were treated equally. The pairwise distances were exported 
from PAUP and were used for the analysis of divergence. To test 
the level of clade support, we conducted 500 replicates of boot¬ 
strapping (Felsenstein 1985) and decay analysis (Bremer 1988; 
Donoghue et al. 1992) using the PAUP program. The aligned 
sequence matrices are available from the first author. The matK 
sequences ot the four species of Liquidambar have been submit¬ 
ted to the GenBank and their accession numbers are AFO 15649 
through AFO 15652. 

Studies have not been done previously to estimate substitution 
rates in the matK gene. Therefore, we estimated a hypothetical 
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Table 2. Liquidambar matK. gene sequence divergence (%). 


Species 

i 

2 

3 

4 

L . acalycina 


0.1 

0.5 

0.9 

L. formosana 



0.7 

1.0 

L. orientalis 

L . styraciflua 




0.5 


rate by using our unpublished data on the divergence between 
two closely related genera of the Hamamelidaceae— Dicoryphe 
Du Petit-Thours (endemic to Madagascar) and Trichocladus Pers. 
(endemic to eastern and southern Africa; Endress 1989)—and an 
estimated time of 100 mya, after which Africa and Madagascar 
were geologically stabilized and direct migration between them 
was probably not possible (Harland et al. 1990; Raven and Ax¬ 
elrod 1974; Schuster 1976). The resulting estimated substitution 
rate is 5.5 X 10" 11 base per site per year for the matK gene 
(assuming a constant substitution rate). 

RESULTS 

Sequence length and divergence. Sequences of the matK 
gene in Liquidambar species were consistently 1512 bases long. 
Sequence divergences ranged from 0.1 to 1.0% (Table 2). In the 
aligned sequences there were 90 variable sites, seven of which 
were found to be phylogenetically informative. 

Phylogenetic relationships. A single shortest phylogenetic 
tree of 92 steps was generated and the consistency index was 1.0 
(Figure 2). The tree consisted of two clades, one ot which in¬ 
cluded the east Asian species Liquidambar acalycina and L. for¬ 
mosana, while the other clade contained L. orietitalis ot w’estern 
Asia and L. styraciflua of North America. Bootstrap percentages 
and decay indices were 99%/5 and 86%/2 tor the two clades, 

respectively. 

Time of divergence. Based on the substitution rate of matK, 
5.5 X 10“' 1 base per site per year (see above), the divergence 
times for Liquidambar species were 9 mya tor L. acalycina and 
L. formosana, 45 mya for L. acalycina and L. orientalis, 81 mya 
for L. acalycina and L. styraciflua, 90 mya for L. formosana and 
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Liquidambar acalycina 


L. formosana 


L. orientalis 


L. styraciflua 
Mytilaria laosensis 

Figure 2. A single shortest phylogenetic tree of 92 steps of Liquidambar. 
Numbers above branches are bootstrap percentages and decay indices, and 
below branches are branch lengths. Consistency index is 1.0. 



L. styraciflua, 63 mya for L. formosana and L. orientalis, and 45 
mya for L. orientalis and L. styraciflua. 

DISCUSSION 

How to estimate genetic differences among taxa has long been 
one of the critical questions about disjuncts in general. Unfortu¬ 
nately, it is very difficult to use morphological, phenological, or 
even cytological characters to quantitatively assess genetic dif¬ 
ference among the disjunct taxa (Lee et al. 1996; Oginuma and 
Tobe 1991). In contrast, molecular data, especially DNA sequenc¬ 
es, have evident advantages over morphological characters in this 
respect (Crawford et al. 1992). Active research has been carried 
out in the last several years, using different sources of molecular 
data to quantify genetic divergence for a variety of disjunct taxa 
(Hoey and Parks 1991, 1994; Lee et al. 1996; Parks and Wendel 

1990; Sang et al. 1994, 1995; Wen and Zimmer 1996; Xiang et 
al. 1994). 

Recent discoveries of fossil flowers attributable to the subtribe 
Loropetalinae (Hamamelidoideae) have extended the fossil record 
of the Hamamelidaceae back to the Campanian (-70-75 mya) 
of the Late Cretaceous (Magallon-Puebla et al. 1996), while flow¬ 
ers of hamamelidaceous affinity extend back to the Santonian 

(~80 mya; Endress and Friis 1991) and the Turonian (~90 mya; 
Crepet et al. 1992). 
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The fossil record of Liquidambar and other altingioid plant 
remains is abundant in the Tertiary, dating back to at least the 
Tertiary-Late Cretaceous, and possibly to the mid-Cretaceous. 
Fossil altingioid fruiting inflorescences of the extinct genus Stein- 
hauera Presl. are known from the Late Cretaceous-Tertiary 
boundary (Maastrichtian/Danian) to the Eocene of Europe (Friis 
and Crane 1989; Mai 1968) and have been considered close to 
either Liquidambar or Altingia (Ferguson 1989; Mai 1968; Tiff- 
ney 1986). Liquidambar pollen was present in the Paleocene 
(~55-65 mya) of North America and northeastern China (Li et 
al. 1995; Muller 1981; Taylor 1990; Tiffney 1986; Wang 1992), 
and Late Cretaceous to Paleocene of China (Sun 1979). Hedlund 
(1966) reports, with reservation, the possible occurrence of Liq- 
uidambarpollenites in the Cenomanian of Oklahoma (~95—100 
mya). 

By comparison, the estimated divergence times (~45-90 mya) 
for Liquidambar species from matK data, using the hypothetical 
substitution rate described above, were fairly close to the tossil 
record, indicating an ancient origin and separation of Liquidam¬ 
bar populations in disjunct areas of eastern and western Asia, and 

North America. 

As noted in the phylogenetic tree (Figure 2), the two clades 
are well supported by bootstrap percentages (99% and 86%, re¬ 
spectively). One clade consists of the two east Asian species, 
Liquidambar acalycina and L. formosana, and the other includes 
L. orientalis of Turkey and L. styracifiua from North America. 
Thus, the matK. phylogeny does not support the recognition ot 
the two sections proposed by Harms (1930) and Chang (1979) 
based on the presence or absence ot setae, and suggests that the 
presence of setae is an autapomorphy. This pattern, especially the 
close relationship of the west Asian species L. orientalis and 
North American L. styracifiua, has been suggested by an allo- 
zyme study (Hoey and Parks 1991). This relationship is also sup¬ 
ported by leaf morphology (Bogle 1986). 

The divergence times estimated trom allozyme data by Hoey 
and Parks (1991). when Nei's (1987) formula was used, were 7 
and 10 mya for Liquidambar styracifiua—L. orientalis, and L. 
styracifiua—L. formosana , respectively. However, when Sarich s 
(1977) and Thorpe’s (1982) formulas were adopted, the times 
tor L. styracifiua—L. orientalis were 13 and 16 mya, respectively. 
The estimated divergence times using allozyme data are rather 
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recent compared to those from matK data and the fossil record 
(~45—100 mya). Although calibrating the substitution rates is 
still not certain for either allozyme or matK. data sources, allo- 
zyme analysis tends to underestimate divergences, especially 
when amino acid substitutions not affecting electrophoretic mo¬ 
bility go undetected (Crawford et al. 1992). We believe that a 
comprehensive study of pairs of disjunct taxa whose times of 
separation have been known more or less precisely, using both 
nuclear and chloroplast DNA sequences, will provide an invalu¬ 
able basis for evaluating substitution rates and divergence times 
of many other taxa. 

The matK gene has been widely used in resolving relationships 
of angiosperms at generic or higher levels (Johnson and Soltis 
1995; Soltis et al. 1996; Steele and Vilgalys 1994), but this study 
suggests that it may be informative also in studying long sepa¬ 
rated species within a genus even though the substitution rate is 
low. Nonetheless, due to the low number of informative sites in 
the matK gene, it seems to be important and interesting to further 
pursue the phylogenetic relationships of Liquidambar species us¬ 
ing a fast evolving region such as the internal transcribed spacers 
ot nuclear ribosomal DNA. 

In summary, the matK phylogeny supports the close relation¬ 
ship of west Asian Liquidambar orientalis and North American 
L. styraciflua, but does not agree with the division of the genus 
into two sections. The divergence of disjunct Liquidambar spe¬ 
cies was probably at least as early as 45 mya. 
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